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Abstract 

We have examined the effect of an oscillatory rotation of a polarized dielectric 
particle. The rotational motion leads to a re-distribution of the polarization 
charge on the surface of the particle. We show that the time averaged steady- 
state dipole moment is along the field direction, but its magnitude is reduced 
by a factor which depends on the angular velocity of rotation. As a result, 
the rotational motion of the particle reduces the electrorheological effect. We 
further assume that the relaxation of polarized charge is arised from a finite 
conductivity of the particle or host medium. We calculate the relaxation 
time based on the Maxwell- Wagner theory, suitably generalized to include 
the rotational motion. Analytic expressions for the reduction factor and the 
relaxation time are given and their dependence on the angular velocity of 
rotation will be discussed. 

PACS Numbers: 83.80.Gv, 77.84.Lf, 77.22. Gm, 41.20.Cv 

Typeset using REVTgX 



1 



I. INTRODUCTION 



The prediction of the strength of the electrorheological (ER) effect is still a main concern 
in theoretical investigation of ER fluids [1—5]. An ER fluid is a suspension of highly polarized 
particles in an insulating host. The ER effect originates from the induced interaction between 
the polarized particles in an ER fluid. Upon the application of an intense electric field, the 
particles in ER fluid aggregate into chains and then aggregate into columns in a short 
response time [3,4]. The rapid field-induced transition between the fluid and solid phase 
makes this material important both for wide industrial applications and for experimental 
and theoretical investigation. 

In deriving the induced interactions between particles, existing theories assume that the 
particles are at rest [6-10]. In a realistic situation, the fluid flow exerts force and torque 
on the particles, setting the particles in both translational and rotational motions. For 
instance, the shear flow in an ER suspension exerts a toque on the particles, which leads 
to a rotational motion of the particles about their centers [11]. Recent experiments gave 
evidences that the induced forces between the rotating particles can be different from the 
values predicted by existing theories [12]. 

To gain some insight into the phenomenon, we have recently formulated a theoretical 
model, which describes the relaxation of polarized charge on the surface of a uniformly 
rotating particle fl3|]. We showed that the rotational motion of the particles reduces the 
induced forces between the particles. We called the reduction of interparticle forces due 
to the rotational motion of the particles the dynamic ER effects ||13|| . In this work, we 
extend the consideration to an arbitrary rotational motion. In particular, we will obtain the 
steady-state dipole moment of a rotating sphere under a sinusoidal oscillatory shear motion. 
We further assume that the relaxation of polarized charge is due to a finite conductivity of 
the particle or host medium. We will derive an analytic expression for the relaxation time. 
The dependence of the reduction factor and the relaxation time on the angular velocity of 
rotation will also be calculated. 
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II. STEADY-STATE DIPOLE MOMENT 



Consider a dielectric sphere under the influence of an electric field E = E z; its induced 
dipole moment is given by: po = PqZ. Assume that it is under a rotational motion of angular 
velocity Co = —uoy. For a rotating dielectric sphere in an electric field, the rotational motion 
leads to a displacement of its polarized charges on the surface of sphere. As a result, there 
is a change of the dipole moment, described by Co x p. The surface charges also suffer from 
relaxation of various kinds, and the rate of change of the dipole moment is described by 
— (p — Po)/t, where r is a relaxation time. The two effects have to be balanced against 
each other, resulting in a steady state dipole moment p, which deviates from the equilibrium 
dipole moment po- Let the resultant dipole moment be p = p x x + p y y + p z z. The rate of 
change of the dipole moment is given by: 

^=Coxp--(p-p ), (1) 
at t 

where the first term on the right hand side is due to the rotational motion and the second 
term is due to a relaxation process, in which the relaxation time r is determined by the 
details of the relaxation process. In component form, the differential equation reads: 

Px Py (Pz ~ Po) 

Px = Wp z , Py = -, Pz = LOp x . 

T T r 

The equation for p y can be readily integrated to yield p y = p y0 e~ J ^\ where J{t) = / ' dt/r. 
Since r (can be time-dependent) is real and positive, p y vanishes as t goes to infinity. To 
solve the equations for p x and p z , we use the complex notation: let p = p x + ip z and po = ipo, 
p can be found by solving the following differential equation: 

dp ( . 1\ ~ . Po / s 

With the initial condition p = po when t = 0, Eq.(|2]) admits a standard solution: 

pe 1 - po = po J ^dt, I = J {^--iuj \ dt, (3) 
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where I is the integration factor. For a uniform rotational motion, uo = ujq is a constant, 
I — - — iuot, Eq.(BD can be solved: 



Po 

P 



(l - iwoTe- t(1_<WoT)/T ) . 



1 — ZU>oT 

As i goes to infinity, we obtain the steady-state solution for a uniform rotation: 

P = i : • 4 

1 — ZUJqT 

In general, analytic solution of Eq.([3|) can be found only for a few simple cases and the 
integral must be evaluated numerically. We concentrate on the steady state solution at a 
sufficiently long time and Eq.(|^) can indeed be solved exactly. By using the L'Hopital's rule, 
we find: 

p = p hm — -r = ; — , (5) 

t^oo Te l J 1 — %{jjT 

where / denotes the time derivative of /. We have assumed that r is real and positive but uj 
can be an arbitrary function of time. Eq.fl^) is the general result for arbitrary rotational mo- 
tion, being of the same form as Eq.(^j). However, the transient solution has to be calculated 
numerically. 

For a dielectric sphere undergoing a simple harmonic oscillation, 6(t) = 6 sin kt, the 
angular velocity is given by u(t) = 9 = 9 k cos kt. From Eq.(^), the steady state dipole 
moment is: 

Po i — 9 kr cos kt , , 

P = 1 - i9 kr cos kt =Po l + 9\k 2 T 2 cos 2 kt ' ^ 

The time dependence of the dipole moment is still periodic. Note that although the sphere 
is undergoing a simple harmonic oscillation, the dipole moment does not exhibit a simple 
harmonic motion. If r is independent of time, we can calculate the time average of the 
dipole moment: 

M = and <gi>= , 1 . (7) 
Po Po Ji + 9lk 2 r 2 



As a result, the motion of particles reduces the ER effect. We define the reduction factor R 
as: 

R=^A = , 1 . (8) 
The reduction is even more significant at high frequencies, R ~ l/d^kr. 



III. CALCULATION OF RELAXATION TIME 

So far, our proposed relaxation time has no explicit expression. If the relaxation process 
is originated from a finite conductivity of the particle or host medium, then we can calculate 



the relaxation time based on the Maxwell- Wagner theory of leaky dielectrics |14j]. For a 
(nonrotating) spherical inclusion embedded in a host medium, the expression is 

/ e l + 2e m \ 

T = e o[ z , , > (9) 



vOi + 2a n 

where ei,e m (oi,<7 m ) denote the dielectric constant (conductivity) of the sphere and host 
medium respectively, eo is the permittivity of free space. For typical values of the permittiv- 
ities and conductivities of common ER fluids, the relaxation time ranges from microseconds 
to milliseconds, and the dynamic ER effect can be observed in experiments. 

In order to account for the impact of a rotational motion on the relaxation time, we first 
replace €\ in Eq.(|9]) by e\ — 1 + where \i is the susceptibility of the sphere. We already 
showed that the dipole moment is reduced by a factor R. If we assume that the polarization 
is uniform throughout the sphere, which can be achieved when the oscillating frequency is 
high, we may write e.\ = l+Rxi in Eq.flj^). Physically it means that the effective polarization 
of the sphere is reduced as a result of the rotational motion, leading to a reduction of the 
effective dielectric constant of the sphere. After some simplifications, we obtain 

r = r oc + . T0 ' To ° , (10) 



where 



, 1 + 2e m \ ( ei + 2e m 

Too = eo — — and r = e 



ai + 2cr m / \ai + 2a. 



' in 



It can be shown that t = tq for k$o = and r — > Too for fc^o — > oo. Eq.([T0|) is a self-consistent 
equation for r and we can calculate the relaxation time self-consistently. 

IV. NUMERICAL RESULTS 

To examine the dependence of the reduction factor R on the angular velocity k0Q, we plot 
R vs. k9o in Fig.|I]for several different values of To. Without loss of generality (i.e., in terms of 
some unit relaxation time), we choose = 1 and r = 2, 4 and 8 respectively. The reduction 
factor decreases rapidly with the increase of k8 , which means that the dipole moment is 
greatly reduced when both the oscillation frequency k and the oscillation amplitude 9q 
become large. 

Next we see how the relaxation time depends on k8 . From Eq.flTOD, r is bounded between 
Tq and Tqo. The lower bound is reached when k8o tends to infinity, which is achieved at 
large frequency. For a larger value of To, the relaxation time decreases more rapidly with 
k6 . The condition of high oscillation frequency reads: 

/c>— . (11) 

The time evolution of the dipole moment is worth studying. In Fig.0, we plot the steady- 
state solution of the perpendicular component p x /po and parallel component p z /po against 
time. We set tq = 2 and r ro = 1, r is then calculated from Eq.([IT]). We first set k — 1 and 
vary the oscillation amplitude 6q. In the left panel of Fig|| we plot p x /po and p z /po against 
time for O — 7r /4, do = it/2 and 6$ = n respectively. For each value of #o> the magnitude 
of the perpendicular component {p x /po) has a maximum value of 0.5 and it has a local 
minimum at t = n/k. When the oscillation amplitude increases, the local minimum value 
decreases, showing a large variation of p x /po- Next we concentrate on the parallel component 
(Pz/po)- When the oscillation amplitude increases, p z /po reduces in general, although the 
maximum value is always equal to unity. This results are expected from Eq.(|7|): 

M = and 1 



Po Po Jl + 6lk 2 T 2 



Hence, on the average, (p z )/po must decrease when we increase the oscillation amplitude 9q. 

Now we examine the case for a constant amplitude 8q but varying frequency k. This is 
a realistic case since in experiment we can hardly increase the amplitude but we can easily 
change the frequency. In the right panel of Fig.jl], we choose 8q = 7r/4 and k = 1 and 3. 
The results for constant k and constant 6q show similar time dependence. However, p x and 
p z show a larger variation in their magnitudes, if we increases the value of k. It should be 
remarked that we have assumed that the oscillation frequency is high so that the relaxation 
time is nearly constant during the motion. 

V. DISCUSSION AND CONCLUSION 

Here a few comments on our results are in order. As our steady-state solution is general, 
one can extend the calculations to an arbitrary rotational motion. We have shown that the 
motion of particles reduces the strength of the dipole moment. It is natural to further cal- 
culate the interparticle force between two rotating spheres. We expect that the interparticle 
force will be reduced substantially because the force between parallel dipoles changes from 
attractive to repulsive when their orientation varies from the transverse to longitudinal field 
case. 

So far, our derivation of relaxation time is based on the mean field theory. We may extend 
the Maxwell- Wagner theory to the polarization relaxation of oscillating particles. In this 
case, we should add a term ppv to the polarization current density, where pp is the polarized 
charge density and v = uj x r is the rotating velocity. However, it is not possible to convert 
the extra term into a dielectric constant and the generalization becomes more complicated 
due to the nonuniform polarized charge density inside the rotating spherical inclusions. We 
are currently examining the solution of the more complicated boundary-value problem. 

The flow in ER fluid may be nonsteady in usual operation situations. But the prevailing 
situation in theory and simulation of ER fluids is to use formulas derived with respect to 
a steady flow. To remedy this drawback, we will endeavor to develop calculation method 
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for suspension hydrodynamics, and use it to study the interaction between particles and 
the oscillating fluid, and derive formulas of the force and torque exerted on particles for a 
suspension [p~5[ |. 

In this work, hydrodynamic (HD) interaction effects have not been considered. However, 
electrorheological fluids are locally very concentrated suspensions and in considering dynamic 
effects, it seems that HD effects can be strong. And this is a future problem. 

In conclusion, we have investigated the problem of how the dipole moment of a dielectric 
sphere varies with time for an arbitrary rotational motion. We have developed a formalism 
for the rotational motion of the sphere and derived the relaxation time by using the mean 
field theory. We have shown that the time averaged steady-state dipole moment is along the 
field direction, but its magnitude is reduced by a factor which depends on the frequency of 
oscillation. As a result, the motion of particles reduces the ER effect. We further calculate 
the relaxation time based on the Maxwell- Wagner theory. The dependence of the reduction 
factor and the relaxation time on the angular velocity of rotation has also been discussed. 
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FIGURES 

FIG. 1. The reduction factor R (left panel) and the relaxation time (right panel) r for = 1. 
The reduction factor decreases drastically for increasing k0Q. The relaxation time reaches its 
minimum value at large kOo- 

FIG. 2. The reduced dipole moment p x /po and p z /Po plotted as a function of time for various 
frequency dependent r (tq = 2, = 1). 
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